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periodic and depended strongly on the amplitude and shape of the surface 
perturbation. In addition to the normal waves produced by wind and shifting 
tides, man-made waves were induced (the bow wave from a small power boat). 


Frequencies generated at the transmitter ranged from 1 to 10 kc. The 





water waves produced by the wind attained amplitudes (crest to trough) approach- | 
at the receiver due to the perturbation displayed several prominent features, 
viz. (1) the shape of the small-amplitude fluctuations matched closely those of 
the corresponding surface waves; (2) the fluctuations of larger amplitude 
exhibited, in addition, harmonics of the frequencies present in the surface 
waves; (3) the average values of both phase and amplitude fluctuations increased | 
with the transmitted acoustic frequency. : 

Figure 1, drawn from Ref. 1, is a plot of average per cent fluctuation in 
signal amplitude vs. wave amplitude. For the lower frequencies (1-4 kc) the 
Signal fluctuation increases with increasing wave amplitude, nearly linearly. 
In general, the slope of each curve becomes increasingly steeper as the acous tial 
frequency increases, and the linear relationship that held for increasing wave 
amplitudes at the lower frequencies is no longer applicable. The results at ; 
5.S kc, which appear to be inconsistent with the other data, were assumed in | 
Ref. 1 to be attributable to variations in water depth caused by tides. 

It was further noted in conjunction with Fig. 1 that fluctuations in 


excess of 30% were seldom observed, regardless of acoustic frequency. 
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In addition, as surface-wave height and frequency were increased, a point 
was reached above which the fluctuation spectrum recorded at the receiver was : 
considerably broader than the surface wave spectrum and showed well defined 


second and higher harmonics not present in the surface wave. In both cases 
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Scrimger's measurements indicated that the significant parameter on which the 


fluctuations depend is the ratio of surface-wave height A to acoustical 


= 


wavelength A (A/A). The value A/A = 18 was found to provide a convenient 
dividing line between those fluctuation spectra containing appreciable second- 
and higher-order harmonic components (A/A < 18) and those spectra containing 
only the first harmonic (A/A > 18). 

In an attempt to investigate the variation in the received signal 
amplitude as a function of surface perturbation position, wave packets were pro- 
duced anc the signal fluctuations were examined as the wave disturbance was 
provagated srom the receiver toward the transmitter or, alternately, in the 
os,osite cirection. In practice the wave packets were created utilizing the 
Dow wave from a small boat on an otherwise calm surface. In order to distin- 
guish the bow wave from motor noise and waves produced at the stern of the 
vessel, the boat was pointed toward either the transmitter or receiver so as to 
avoid cutting the imaginary line connecting the two instruments, and the motor 
Was stopped several seconds before the wave train reached the tripod. 

The experiment was carried out initially using a single transmitted c.w. 
frequency. Later, in order to provide reliable results on the effect of 
changes in acoustic frequency, two c.w. sound signals were transmitted 
Simultaneously. 

Wave trains were propagated from receiver to transmitter and subsequently 
rrom transmitter to receiver for five frequencies ranging from 0.7 to 2.0 kc. 
The spectra recorded at the receiver indicated that fluctuations occurred only 
when the wave packet was in the immediate vicinity of either the source or 


receiver. As higher acoustic frequencies (up to 10 kc) were introduced, fluc- 


tuations in tne received signal appeared when the surface wave was located at 
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greater distances from the source and receiver, but the maximum disturbances 
continued to occur when the perturbation was directly over source or receiver. 

A frequency was eventually reached (approximately 10 kc) which produced fluctua- 
tions at the receiver independent of the position of the wave packet. In all 
cases, the presence of harmonics in the signal fluctuations depended on \/A 


as mentioned previously. 
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A summary of the conclusions drawn from Scrimger's paper follows: 

1) The average fluctuation in received signal intensity increases with 
increasing frequency as well as the increasing wave height, up to a maximum 
rliuctuation value of about 30%. 

2) For large values of A/A (i.e. acoustic wavelength to wave height), 
the shape of the fluctuating acoustic signal closely resembles the surface-wave 
shape, whereas small values of A/A go with received fluctuations that are 
broader and that contain more harmonics than the surface wave. A convenient 
dividing line for \/A separating the two cases was found to be approximately 
es 


3) Propagating a wave packet over an otherwise calm surface, between 


mainly while the perturbation was near the source or receiver. However, with 
increasing acoustic frequency, the fluctuations were also observed for wave 


packet positions further removed from the source and receiver. 


source and receiver, showed fluctuations in low-frequency signals to occur 
{ 








As a first step in providing a theoretical explanation of Scrimger's 
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SECTION II -- Development of the Mathematical Model | 
| 
| 
results, the omnidirectional cw transmitter used as a sound source throughout 


the actual experiment can be replaced by a simple point source located a dis- 





tance h below the water's surface. In an unbounded medium, this source spreais 


l 

t 
energy uniformly over a complete spherical surface at each radius "r" from the : 
origin. The acoustic pressure P is given by the expression | 
( 
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iCwt-kr) 


P = [ip ckQ./(unn) Je (2.1) 
wnere 
| 
! c = sound speed in the water 
| v = frequency = ck/2n 
W = acoustic angular frequency =217Vv 
k = wave number = w/c 
Daun equilibrium density of water 
r = distance from source to desired point in the sound field 
: oa = strength of the simple source. 
In addition, the problems encountered in dealing with the air-water 


effect between two sources: the original source, given by Eq. (2.1), located 


a distance h below the reflecting surface, and an image source of equal 





strength reflected in the air-water interface and located a distance h above 


the surface. The acoustic pressure for the image source is given by 





interface can be simplified by considering the acoustic field as an interference} 


: a E. Kinsler, A. R. Frey, "Fundamentals of Acoustics," (John Wiley & Sons, 
inc., New York, 1962, 2nd ed). 
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p = ([- ip va /(ar) Je (2.2) 
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where r is now the distance from image to receiver and the minus sign, which 





indicates a 180° phase shift at the surface, satisfies the boundary condition 


that the total pressure given by adding Eq. (2.1) and (2.2), must vanish at the 


= 


surface. 


bottom have been completely ignored, we have a dipole expression for the pres- 


sure amplitude at any point in the sound field (see Fig. 2); 


-ikr, -ikr, 
Pa = [ip va /2][le ye - (e /v,], (272) 
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Combining Eq. (2.1) and (2.2), and recognizing that reflections from the 
distance from point source to field point. 


2; 
Si 


distance from image source to field point. 
: 
It can readily be shown for spherical waves that the radial pressure 


gradient <dp/sr is related to the radial particle acceleration by 
- dp/or = p3°E/ate 


where the symbol §€ represents radial particle displacement. 
Orienting the coordinate system as indicated in Fig. 2 with the origin . 

located at the point 0, and considering only the y-component of particle 

motion (i.e., perpendicular to the air-water interface), we obtain the 


particular case 
- = 2 2: | | 
dp/dy p9 E/ot er en 


where y is the particle acceleration in the y direction. 
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Figure 3 
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It is evident from Fig.) 233 jeu and r can be expressed as follows, | 
, 


in terms of x, y, and ih: 


[  yj217!? (2.5) 


Pio = [x? + (h + y) 


Substituting these expressions in Eq. (2.3), and using Eq. (2.4), yields | 


oe -ikr, -ikr 
[Ce 
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ay = a/dy{(ivg /2)e - (e *Y/ry]} 3 Cer 


-ikr -ikr 


‘Ot Tik(h-y)e  *I/r2 + [(h-y)e — *]/e3 


a = (-ivQ /2)e 
y O 
-iky, -ikr, 
+ [ik(h+y)e liar + [(hty)e coe (25 
betting y > 0, 1.e., considering only those partacies locatedyacs tne 


air-water interface, we obtain from Eq. (2.7) 


a lwt -iks 2 . -lks 3 , 
ae = vhQ.e BL Gke )/s (ie )/s*], eae 


where s = ee the value of rv, and vr, when y = 0. 

Thus far, the surface has been considered to be plane and free of waves. 
As a first approximation to the surface-wave perturbation, let us use a single 
symmetric pulse moving along the x axis and instantaneously located at x, 
és illustrated in Fig. 2. It is physically obvious that such a distortion in 


the reflecting interface will alter the acoustic field at the receiver. It is 


also obvious that this new effect vanishes as the height of the surface-wave 


assuned to be very small compared with the distance from source to receiver, a 


plausible epproach is to regard the disturbance in the acoustic field as coming 


| 
pulse goes to zero, and, since the horizontal dimensions of the pulse are 
| 
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| 
1 from a dipole, with a vertical axis, located at the pulse. The main task is to 
determine tne strength of this dipole. 

: ae 

7 Letting Pl! vyvepresent the incremental acoustic pressure caused at (x,y) 

i v 
| by the surface-wave ~- which is assumed to act as if it were a dipole -- we can 
woite i 
ie -ikr? ~ikr, . 
I « ~ 

| P! = [Cip_vQ')/2Je” [le \/ei - Ce ia (2,938 

| D oO oO ahs 2 

! ¥ 
( where the perturbation-dipole strength Q‘' depends, in a@ manner to be deter- | 
| 
7 minec, on the source-dipole strength Q, and the source separation d!' for the 










)  senturbation cipole is proportional to the wave height (d' enters into ry 


ana xm, just as d= 2h did in yr, and r,). Henceforth P as given by 


L D? 


Zc. (2.3), will be referred to as the “source dipole," and Pi» as given by 
. (2.9), will be referred to as the "perturbation dipole." 
Figure 3 shows the distances Recied to specify the pressure amplitude 
) occurring at the receiver for an arbitrary position of the perturbation aapones 
If one calls the mass of the surface-wave pulse m, the vertical force 


i: exerted by the field of the source dipole is given by 


i 2,8 mvhQe*"*[ (ke ***)/s2 Ge ny ee Fie, (230 
where B70 is determined from Eq. (2.8). 


| In the preceding paragraphs we have succeeded in specifying the particle 
. 


ecceleration due to the source dipole Pos anywhere on the surface of the 





air-water interface. In addition, through knowledge of the surface-wave mass, 
the force exerted by the original acoustic field upon the perturbation is given 
by Eq. (2.10). We now attempt to express the perturbation-dipole strength Q' 


meecerms otf this force. 
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The fluid volume in the wave pulse is subject to the dynamics of moving 





compressible fluids and therefore certain conservation laws hold: 


Ie, 2ne Equation of sConpinuas, 2 
s0/3t + Ve{pv) = 0, Qa) 


where v is particle velocity. | 
This equation expresses the fact that for a unit volume there is a balance! 
between the change in density and the masses entering and leaving the system 
per unit time. 
A source function can be inserted into the continuity equation to account 


Zor the rate at which mass is introduced. Thus Eq. (2.11) becomes 
dp/ot + Ve(pyv) = 64> (2a) 


where 


oD 
tH 


instantaneous density 


equilibrium density 


OD 
if 


strength density. 


O 
51 


Since the spatial variation of density is much less than the spatial 


variation of the particle velocity, Eq. (2.12) reduces to 





dp/at + Oh Sapna (2.12a) 

II. Conservation of Momentum -- 
! 
po (ay/at) + po (yrVy) + Vp = £. (2.13) # 
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Equation (2.13) is simply Euler's equation of motion for fluid flow, wit 





Mumene presstme, and f£ the external force exerted on the fluid. The mass of 
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Fluid beinzs considered is that in the su.face-wave pulse, which is located 
essentially at y = 0. Equation (2.7) shows that for y = 0 the particle 

f 
acceleration ve is indepencent Ore ees ee ee which is mere) atte | 


time integral of oe is also indepencent of y. Hence v has at most only 


Par ete 


the particle velocity v is a only -- i.e., particles move only normal to 


} 
x- and z-components. But we are assuming that at the pressure-release surface 
@eeour-ace. It follows that in the perturbation pulse v-Vv vanishes jean 3 


ey 
)O 
co~ 
tO 
}— 


{ 
3) seduces to | 
po (av/ot) + Vp = f. (2.13a) | 

O fv ! 


Taking 9/dt throughout Eq. (2.12a), operating with (-V-) throughout 


aq. (2.134), and adding the results gives 
eo/at? - V2p = p (aq/at) - Vg. (2.15) 


We assume that the fluid is homogeneous and nondissipative, so that any 
Change 4p in the density is proportional to an accompanying change Ap in 


the pressure: 


Ap = Ap/e?, (2.15) 


where c is the speed of sound in the fluid. It follows that 949/a¢t can 


Berrepieced by Cao ye and Eq. (2.14) becomes 


c72(92p/at?) - V2p = p6aaq/at) - V-E. (2.14a) Ff 





Since the single pulse representing the surface perturbation is assumed 


to be very small compared to the area of the calm surface, the external force 
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exerted on tie wave by the source cipole can be approximated by a deita- 
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mese > «eS cCOomCcnen Ts x, Y, Zge anc tne center of the wave pulse iS vageee HP 
“ny = ““"O 
che stremsth density «a in =q. (2.122) can be written in the desired | 
Elpoie form, as | 
q = Qc VLéla-9 dds (2.10 @ 


where Qd is the dipole strength of the perturbation dipole (Eq. (2.9)) and 
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d is the vector separation of the two point sources in this dipole. Then 


Eq. (2.i4a) becomes 
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Veto (30/dt)d6( rn gee OC eee (2.18) 


The last step is allowed because (V+) does not operate on the constant Po» 
nor upor. Q and d, which are assumed to be independent of x, y, z (i.e., 
the perturbation wave does not change shape as it travels, and therefore the 
essociated perturbation dipole does not change shape, either). Also, (V+) andi 
9/3t can be assumed to commute. | | 
Equation (2.18) has been developed for the general case in which both a 


generating source q and an external force £ are present. In the true physi- 


cal situction of this problem, there is no source putting fluid into the wave, 


Gir 


but there is an external force fF produced by the original source dipole. Then 


the Sight-nend side of Ea. (2.18) reduces to | 
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pulse, woulc be unchanged if we postulate a nonphysical source whose effect in 
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Therefore Eq. (2.18), which governs the acoustical behavior of the surface-wave 
mee. 2.18) is identical with that of the ectual external force f: 
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Evidently d, the vector separation in the nonphysical dipole, must be parallel 
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Cwm y-components. From Eq. (2.16), F equals Pies end to be consistent 


With the notation in Eq. (2.9) we can write 
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Q=Qte". (2.22) 
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Gee 2ing 7, from Eq. (2,10), we finally arrive at an expression for Pa of 
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wawre the unperturbed (and major) part P. comes from Eq. (2.3) and the small 
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2.24). The wave height appeers implicitly in P 
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Meagicications ane introduced in order to acaot the model to the physical 


In the actual experiment verformed by Scrimger, the per cent signal 
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\ feet oations were obtained directiy from magnetic tepe and strip charts usedaro 


record the variations in pressure emplitude. Straight lines were drawn through 


' 
the upper and lower average level oF the peaks in the fluctuating signal to 
the mean signal level I by the expression AI/2 x 100/I to produce the per 
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Fluctuation. We shall hereafter omit "100," using the relative fluctuation 





Although, as noted above, the original article expresses the per cent 
Tluctuation in terms of pressure amplitude, it was found more convenient to 
| consider this percentage as a function of signal intensity. For small values, 


my 2AP/P, 
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ressure amplitude at the receiver due to the source dipole is given 


give the quantity AI for the over-all fluctuation. This was then related to 
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the second bracket, the binomial theorem can be employed and Eq. (2.31) reduces 
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As in the cevelopment leading to Eq. (2.29), replace r‘, a in the 
-., 


Bemo.inatcros Sy 2I5 and in the exponentials by Pio chd'/2n, and Yi othd' /2r, 


respectively. Carrying out the same procedures used to obtain Eq. (2.29) on 


jeewcemainings exponertial factors, we find for PA 


: ~ik(v! +s) 
Oe ie ae 7etgees sin(khd'/2r 
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fr Soth Eq. (2.27) and Eq. (2.33), the fact that the second term ie mucus 

smailer than the first is not sufficient reason to employ the binomial expansion 
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feeeeral! r's are multiplied by k to form the exponents ¢é . The 


important test is that the next term in the binomial expansion — omitted here — 
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Defining the intensity I = CP P, where * indicates "complex 


conjugate," C is a dimensional constant of the same value for Py and Pas 
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and 


P = total pressure amplitude = P, + Py, 
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Naving neglected the tern ee a 


(2.36) 


In evaluating AI/I, we note that the factor p,vae"* is common to 
bota  P) (i.@., Eqae (2529) ad Ph (i.e., Eq. (2.35)) and so will cancel 
Memea. (2.36). Applying Eqs. (2.29) and (2.35) to Eq. (2.36) and oncemacaud 

aa Gr 0 ¢.. see ; 
uUtizising the fact that e -e = 21 Sin Gee, we obtain 
Vkhr Sin(khd '/2n" =) 
a= ga NS See | sinlk(r}, +S - Pi5)4 OST) 
d'ar!,s* sin (kh*/r.,) Ss 
12 ji 
where V=m/p. is the volume of the wave pulse. 
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SSCTION iii ~~ Evaluation of the Proposed Model 
An examination of Eqs. (2.37) and (2.24) will help us to compare the 
benavior of these expressions with the results obtained by Scrimger. 


The intensity fluctuation is given in Eq. (2.37) as 


+ 


Vkhr Sintnnd 2a" 
hie ee sinlk(r) 5+ 


t 24 ce Zin 
Tr) 5S d sin (kh [24 5) 
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Equation (2.37) would be much simpler if each, or any, of the three 
trigonometric functions couid be replaced by their respective arguments. We 


consider first sin(khd'/2r),). When the wave packet is directly over the 


receiver, Pi 5 takes on a@ minimum value which is practically equal to the 
esa , 1 : : 
water dept. h, and the argument of the sine becomes = kd . Figure 1, which 


is crawn from Scrimger's paper, shows that the maximum product of wave number 


K <vimes wave height is 1.1. Since d' ovresumably approximates the wave 


Vy) 


a@eezht, this sine function can be replaced by its argument. Thus Eq. (2a 
becomes ae 
Vk«“h“r 
AI/I = - [= Sasa eee sinlk(r! See eal (2.1) 
Di es See a eee 12 * me 
12 12 
and since d' makes only a negligible contribution to Pio in the last sine 


term of Eq. (3.1), and appears nowhere else, it follows that, for kd' < 1.1, 


AT/i eave (3.2) 


"Several points corresponding to 5.9 and 10.9 ke exceed this product. However, 
the points associated with each of these frequencies are inconsistent with 
the remainder of Fig. l. 



































Consider next sin(kh7/>,,), Ini WelC hale paeeana Yi5 are constants such 
chat sin tkh*/r,.) becomes approximately sin(k/3). To replace sin(k/3) 
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ereater than 0.7, thus making 


x < 2.1 ft” (frequency 1.7 kc). For this value of k/3 the discrepancy 
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| for =he sine, Eq. (3.2) becomes 
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last remaining sine term, AI/I varies with frequency in a linear manner. As 
) x/% vises to about 7/2, sin(k/3) is nearly independent of frequency, and 
“I/= becomes proportional to w*. When k/3 +7, sin(k/3), as well as Pos 
approach zero (see Fig. 4). The fact that AI/I approaches 
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| Thus at low frequencies (i.e., less than 2ke), and still ignoring the 
Meetaity wher. sin(k/3) approaches a can be attributed to the fact that the 


receiver is in the Lloyd-fringe region where small changes in frequency, or in 
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1 (tides, etc.), will produce large changes in the fluctuations measured at 
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| the receiver. 


Still postponing an examination of sin k (pl 5tS-P)5)> consider the 


125 found in Eq. (3.3). In it, r', can safely be re- 


eZ 


and Fig. 3 shows that s = Coon Therefore we obtain 


| Pypmactude term r 


placed by Yi> 


BS yss- Ce) ae Caney a (3.4) 


| where the expression for r! 


; is obtained from Eq. (2.34). Temporarily chang- 


| ing the coordinate system so that the origin is located midway between source 


W and receiver (i.e., x= 0/2) and letting n= x-r,/2 alters Eq. (3.4) to 
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rance poinc (see Fig. 5), even for values of x to the left of the source and 


@ewene right of the receiver, as seen in Fig? ae 


Taking -dB/dn @in Eom. 5) vyieice 


ab dna 2n? ~ 


ES 
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+ Qnh?, (3.6) 
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anc setting dB/dn = 0 will disclose any maxima or minima in B (there are 
three such points since dB/dn is cubic in n). Further, by taking 4d*B/dn* 
anc investigating the sign of the resultant expression, one can establish 


wnetner & maximum or minimum exists at any given n satisfying dB/dn = 0. 


a) At n = 0 ~- mid-range -- B has a maximum. Substituting the 
proper numerical values for Ys» h and n into Bq (3ee)) produces 
B= (= v2+h?)? % fie ~ 5 x 10°Ft*, ee 


b) At n= (} 2-n2y1/? -- i.e., nearly over the source or receiver 
but displaced slightly toward the center of tne system -- there are minima, and 


Eq. (3.5) becomes approximately 
BR Doe ta \ Oe OC Ete | (sac) 


Thus, still without considering sinlk(e),+s-r,5)1, Eds 4 cmc erence 


and (3.8) show that AI/I is about 50 times smaller when the wave pulse is at 


che mid-range position (n= 0) than when the wave pulse is approximately 


over the source or receiver (see Fig. 5). 
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Again investigating Eq. (3.4), assume that x <<r although not 


i2 


| necesSarily << h (i.e., the wave pulse is in the vicinity of the transmitter).]} 


He@memion (254) becomes 
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BX (x@th? )(r2+h?) & 2o(x?th?), (3,9) 
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since ve ->» h2. Thus the value of B is approximately symmetrical about the | 
transmitter (at «x = 0), and by a similar investigation over the receiver 


(rn, -x << 21). the same approximately symmetrical behavior is observed. 
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The last step is to analyze sinlk(r!,+s-r 
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wnere Ys 5 ence i. 12 have been defined in Section II and gs = (2a 


In order to investigate the behavior of Eq. (3.10) we consider three special 


cases; [|x| << vr, me -x| = jul << Se) 5 |x ae ahs niece D 


5 In Case Dpiveea 
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x becomes small, r! is nearly equal to r and the term hd" /2r5 in | 
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Eq. (3.10) can be written fairly accurately as hat /2r, 
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peur th?) eae or pe/(epthe)] + hd'/2p, + (x2+h2) - h¢/r, Cay 
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where x* appearing in the first brackets of Eq. (3.11) is considered negli- ! 


gible in comparison with 27,x. 


Using the binomial expansion on the term in brackets of Eq. (3.12), and 
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multiplying by the first term, we have 
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| Again using the binomial expansion on oorn im Eqs to ee 
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approximating (rs7+h*) wr. in the small terms that resul tye «cusps 
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= (x oho = n*/(20,) “- hd'/(22,). (3.1l4a) 





Thus, for very small x (i.e., near the source), sintk(r},ts-r,,)] is 


given fairly accurately by 
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sintk{ (x?+h*) x - h*/2(r,) + hd'/(2n,)]}. (3it5 
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The other cases go similarly, after coordinate transformations from x to u 


meet} are applied. Thus, Case 2, for small wu > x—=e.) where a Aas 
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Case 3, for small n = x- =P, (i.e., mid-range), leads to the form 


sin{k[h?/p, + hd'/r, + hd'n/zj Jt. (saa 


Comparison of (3.15) and (3.16) reveals that except for the terms in 
d' they behave symmetrically in their respective areas of influence if account 
is taken of the signs associated with x and u: x > 0 means going from 
source tcwarc receiver, whereas nu < 0 means going from receiver toward 


source. EBotn terms in d' are small (d' is presumably a few inches at most). 
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Because of this similarity between (3.15) and (3,16), it suffices to examine 


f (3.15) only. 
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When x goes to zero, (3.15) becomes 
sin{k{h - (h/2r, )(a-d*) J} ~ sin(kh). (3,18) 


When n-> 0, (3.17) becomes 
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sin{k[a*/p, + hdt/r, + hé'n/ep]} % sin(kh?/r,). (2.19) 


Thus, even at the lowest frequency, 0.7 kc, considered here (correspond- 


ing to k%i), the argument or this sine term changes from about 10 to 0.33, 


, - , at 
Mee oO. -cOut 37 to about O.in 3. 2S - Kee geen oO FeO a Ty: As 


tne argument of the sine practically reverses 
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x gecs Powe oO > =a eo one 
é, The change is proportional to k and therefore to frequency; 


Suememecdc OS Kc the argument or the sine chenges from about = 207 “down te 27 
ema Dack to ebout 207 , as x traverses the range rrom O to ro: Hence 
this sine term, the last factor in Eq. (2.37), can never be replaced by its 
ersumenc. It modulates the overall change in AI/I produced by the variation 


with x, making AI/I become alternately positive and negative 
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In the region around mid-range, Sars ae is so small that it matters 
little what the sine term does, but near x = 0 and x=pr, the effect of 


intk{h-(n/2r, )(h-a") J} is important. 

Since h and yr. are fixed in this investigation, the exact value of | 

* e | 

an ix(pr +s->,,)} at x= 0 or r, depends mainly on k, though to a | 
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Prat. extert on d*. The question arises how rapidly does this term 


chenge near the source and receiver, where ~- because of the small a a 


AI/I can be large? / 
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To obtain an approximate answer, we neglect the small terms sia 2s 


and hd'/2r,, in Eq. (3.15). Regardless of the value of the sine at x = 0, 


it will be sure to have a zero when the argument changes by +1. Hence, to 
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e in x we write the argument of the sine for x = 0 
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27h/ ~ (Qn /, i (n2en2)7*!? j= sae (3526) 


Mee. Kk = 27/A. Solving Eq. (3.20) for x yields 


|! 
i 
i | i | 
x © = ACL /4n) or xn - = AC(1-A/4R), (S624) 
Bee pulls @id minus 7 respectively. 
aus feom Eq. (3.21) ana Ea. (3.15) 5. te the emreme it Can Dev irueres 
(i.e., only at smell x), it is seen that when the wave pulse travels a dis-~ 
cance A, éepproximately but not exactly centered over the transmitter, 
cin (kl (n2n2)-/? nh? on, the hes at 1 ros ( e (2.15) | 
en ( xX“ 4 —KH2 oth L208) as at least two zeros (three, i 
hepoens to be zero exactly at x= 0), and AI/I changes by a full cycle. a | 
Similiar argumert holds over the receiver, because (3.15) and (3.16) behave | 


So nearly identically. In contrast, (3.17) shows that near mid-range the sine 


function scarcely changes with n, and so with x, bur AI/I is small in 


tnet region, as we have noted, | 
AT/=. as discussed depends -- except for its amplitude and part of its 
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[meas well as for small + Ax near x = 0, end x - Moreover, an examination 
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of the conclusions drawn from (3.15) through (3.17) show that kL?) FSS 
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<s nearly symmetric about the mid-range. Since Ps is a constant throughout 
fee problem, Kp) ts) must also possess this symmetry. Therefore we conclude | 
Sie tS! . 
inact e 7 » which represents the phase angle associated with the pres- 
| 
sure amplitude of the perturbation dipole, is aiso approximately symmetrical F 
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Thus if a line is drawn perpendiculer to the air-water interface and 


4. 


located micway between source and receiver (i.e., mid-range), the pressure 
amplitude cetected at the receiver due to a surface wave located to the right 
of this line will produce nearly the same phase, at the receiver, as any wave 


Mowume corresponding position to the left of the reference line. | 
Mme snort , Py has much the seme properties as AI/I. As a test of the 


low-frecuency approximation introduced to obtain Eq. (3.3), some exact calcula- 





tions of | have been made from Eq. (2.24) for several frequencies. The 


results are shown in Fig. 7 (pp 36, 37). 
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Before starting detailedcomparisons of the present theory with Scrimger's 


measurements (Ref. 1), we point out thac related buc somewhat different quanti- 


he two cases. Most of Scringer's conclusions present the time 
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aonituces of measured fluctuations in the received acoustic 
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pressure, converted to percentages or the magnitude of the average pressure 


“memeer, i= Cris thesis it has been more convenient to deal with the vaive 
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AT/I vrocuced when the perturbing wave pulse is at some specified position x. 


In Scopinges's cise, many waves and trains of waves have passed by source and 
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PecSiver for cach plotted value of fiuctuation. Hence we should imagine our 


Weave zulse to passover the whole ranze -- contributing strongly to fluctuations 
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@verace of AI/I for values cf x near 0 (or near oi). The corresponding 


+s~r,_,)] by the average of its magni- 





tude,.i.e., by 2/%, There is also, of course, a factor of two between AI/I 
eue@ @-/P, end of 100 in converting to percentage, but these are constents and 
of no interest unless absolute magnitudes were to be compared. 

The sine factor just mentioned is not always irrelevant. It would, for 
€xample, play an imporcant parc in discussing the frequency of individual fluc- 
tuuvions es a weve passes by. Therefore we prefer not to average over it except 
when this step is dictated by the nature of the experimental data. 

Sauvetion (2.37) proved reasonably successful in reproducing the experimen- 
tal results determined from Fig. 1. An examination of the points making up the 
al pram in the Scrimger peper, upoh which Fig. 1 is besed, shows that for 


Pe@cwencies bellow 4*kce the fluctuations appear to have a Linear dependence on 





Fmeeve neicht. For frequencies at and above 't kc, this Linear dependence becomes 








= a a _ 
ss ee = ~ I mn ——— — 
Care lee il ca Vet TIS Dai a LE enti al Bal Elin Di LF TT ee Mian DR eabealnat Ve bed esti bi ileal ebb eenemchse Deusen Sea Soazsous ene = - = rt —_ ~ © a 
aera bina = ae : fit tal si intiaen Septal De , Seam i i Se gl Dine IE iene Silent a 2 si ie TE Rj 
D 
i 
—— 


questionable, the data points being quite widely scattered. As noted in Section 


I III, AI/l is proportional to the volume V of the water~wave pulse. This || 
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mi) volume is certainly dcepencent upon tne wave height, since increasing tema. 
| 
: ¢ certainly increases th it is necessary +t >e whethe j 
height certainiy inereases the volume, but it is necessary to see whether this 
| GCapencence is linear or not. t 
The puise wave treated in this thesis hes, thus far, been only partly 
{ 
specified. It is quite smail in horizortal dimensions -- e.g., much smaller : 
i 
than m, oF h, and presumably not large compared with an acoustic wavelength; 
co 


related cisole sepération d'. To make comparisons between theory and Scrimger'’ 










measurements, the wave pulse should be taiiored in shape to resemble a piece of 


the sort cf weves he used. Scrimger's rigs. 1-9 are associated. with various 
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also, it hes a@ height proportional to end presumably not much different from the 


trains of meexly pavallel wave crésts, with small heights, presumably extending 


€@ consideresle clevance. In any ore test, the weve train had a fairly definite 





Beeeecicy v , rrom about 1.0 to about 2.5 sec ~ in various cases. Therefore 





nhac a fairly definite wavelength A which cén be computed from 







v g/ (212 ) | (4.1) 






where g¢ is the acceleration due to gravity. The amplitudes of the waves were 





always rather small compared with A, and therefore the waves were close to 





: ae by 
a sinusoicai shape in vertical section: 






u(x} % =A cos(27x/A). (4,2) 





EvyCoCamee™. ies. "CD Dover Publications, New York, 1945, 6th ed.), p. 417. 
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yea of a part of such a wave is easily shown to be directly 
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proportional to A (and to A as well, but this is not relevant at present). | 
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Hence we can define the shape of ovr wave pulse as follows: In the ; 


direction parailel to a wave crest, it will have e short constanc length, not ) 


a function of x,A,A, etc.; in the vertical plane at right angles to a wave 


crest, it will be shaped like pert of a sinusoid. Then the volume V, as 
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used in previots eauetions will be proportional to the wave height, or to f 
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12S been discussed separately -- the relative fluctuation in j; 
B@ecnsity is proportional to kA and so to A/A. Thus the linear dependence | 
Of the fluctuations on wave height observéd by Scrimger at any one frequency ; 
i 

below 4 kc (anc possibly for frequencies above this value) is confirmed by the 
Pescuics chtadinea in this work. | 
Figure 4 is a log-log plot of relative signal fluctuations vs frequency, | 

Trom 0.7 to 6.5 kc, where the ordinates are proportional to the logerithms of 

[AT/I)/V}] from Eq. (3.1). The theoretical curve in Fig. 4 is linear in fre- 


G@uency meer v = 1 kc, increases to ae quadratic dependence in the vicinity of 


is attributable to the fact thet for high frequencies the receiver is still in 
tne Licyd-fringe region, and near 7 ke PD should theometically go (to zeroe 


making AI/i approech infinity. In the neighborhood of 7 kc, small chenges 


v¥ = 2.5 kc, and finally approaches infinity as v #7 ke. This last phenomenon 
i 
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in frequency or in water depth will produce large changes in the fluctuations f 


obtained at the receiver. The experimentally obtained curve of Scrimger is 
i 
also plotted in Fig. 4, for frequencies from 1.1 to 3.3 kc. It shows a nearly | 
Seearasc.— depencence upon frequency as v shifts from 1.1 toward 2.0 kc: es 
increases from 2.0 to 3.3 kc, this dependence becomes more nearly lineer, 
ne everaze siope between 1.1 and 3.3 ke snowing a dependence midway between | 
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linear end quadratic. aus in tae region Setween v = 0.7 end 4 ke the frequency! 
t 
Cepencence of AT/T nvvredicted by our mathematical model is in fairly good agree 


ment with Scrimger's experimental results. For frequencies higher than 4 ke 
this egreement breaks down; nowever, above 4 ke the reasoning used in Section ITi 
to replace Sinha terms by 1tS argument) eiso loses Pes velicia, 
Fisure 5 is a semi-log plot of lao ~ vs distence from the source, 
to show the unmcdulated symmetry present in Eq. (3.3). Figure 5 investigates 
~ 26 PA = fa = iS ‘ Sees earatoe abe i 
(x 95 ) ~ from mid-raenze (i.e., 150 ft to the right of the source) to a point 
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fe. 37). the form of AI/I @s 2 Function 
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between 285 fx and S1l2 ft, to ellow shewing sufficient detail. Again, results 
computed for a similar region about the source produce identical results; the 
efrect of Ec. (2.37) ebout the mid-range is negligible as noted earlier. The 


lighter lines (i.,e., the envelope) represent 
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A varies from 1 to 5 rt for the freauencies under consideration, and one would 
expect AI/I to show a mixture of the two rates of variation. Finally, 1t is 
well to remember here than when considering a continuous wave train the total 


Po becomes a phase-dependent sum of several Py's from individual waves 


located over both the source and receiver, thus further complicating the measure 


maeccerns. 


Figure 7 is a semi-log plot, vs horizontal distance x of the wave pulse 


from the source, of the quantity BRE 


~ Bea as calculated from Eq. (2.24) 
before various approximations were introduced. In the calculations a value of 
1.14 inches was chosen for d', intermediate between the extremes of A given 
maeRer. 1. he absolute values of |P|* are unknown here, since they depend 
on both oe and V; but the correct relative heights are maintained in all 


cases. The Figure serves as a partial check on the approximations, introduced 


earlier, to obtain Py and AI/I. The following points are noteworthy. 


a) "| * behaves approximately the same whether the wave pulse is over source 


Or receiver, but more nearly the same at low acoustic freauencies than at high. 


fb.) et“ is more than TO times as great for x ¥%0 or a than ior 


: 
x *% ay (the approximations of Section III would predict about 50* or 2500). 
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(c) |p! 


D increases for all x with increasing acoustic frequency. (d) To 


some extent there is confirmation of Scrimper's observation that at higher 
frequencies fluctuations are observable even when the surface perturbations are 
near mid-range. The ratio of maximum to minimum in Fig. 7 is practically inde- 
pendent of acoustic frequency, but the minimum -- as well as the maximum -- 
becomes much higher as the frequency increases. Hence fluctuations caused by 
Waves near mid-range might be undetectable at low frequencies but detectable 


atu mign ones. 
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\) In summary our theoretical results are in fair to good agreement with most 
f 
} 
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' of Scrimger's measurements, for the lower acoustic frequencies, as outlined ; 
| | 
| 


; below: 
i 
tf ; ; k 
: 1) When the wave pulse is asnroxivately above the source or receiver, i 
| ! 
SeAL/I will se far larger than when tne bulse 25 Im the mid-range regiene ft 
Hi 
2)  .« Fixed acoustic frecuency, 4£./2 will be proportional ico sia 
i 
Simei soc ffl 
) . 
3) At fixed A, AI/I will rise with frequency at roughly the same rate 
| 
| omic 4 SCrimger Sound (Fig. 4). : 
4) More generally, AI/I is a function of A/d, rather than of A or 
| ; 
; 
| A Beoaratc y. 
‘ 
. f f 
' In tnreée other respects, the present theory is not complete enough to be 


Be usea cor comparison with measurements: 


| 5) «agreement in absolute magnitudes of AI/I is not to be expected, 
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W uMt=l smail weve pulses as considered here are joined into trains of several 


t.Se 
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| Sardilel weaves extending a considerable distance in the z direction (across 


the vange). 
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6) Any effects peculiar to high frequencies -- e.g., the saturation shown, 
mamecrimger's Fig. 3, and breakdown of the rule that only waves nearly over the 
source or receiver are important -- cannot be treated in the present approxima- 
tion for low frequencies, although Fig. 7 partly corroborates the latter 


example. 
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| 7) Finally, it is not evident why harmonic and other distortions should 
appear in AI/I (Serimger's Fig. 9) when A/A rises above the approximate 


value l/lée. “rite possibly our assumption that the wave pulse acts like a 


ee ce eR Pt PEER, I A 


Simole accoustic dipole fails when the wave height is too great, but this is a 


SPisuetise anc rot an explanation. 
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In an attempt to cevelop a rudimentary and workable modeiueapas mma. | 
recreating the experimental conditions, several parameters and boundary condi- | 


tions were ne 


glected or greatly simpiificc. These modifications, and suggested 


refinements of the model, are listed belcn. 
A. Shave of Surface Wave -- In the development of the problem carried out | 
a i 


f| 

in Section Il, and ciseussed in Section III, the surface-wave perturbation was |] 
| 

@eMer bed 2s a single symmetric pulse moving from scurce to receiver. In actueum 


practice, a surface wave, whether natural or produced as the bow-wave disturbance 
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Beem 2c bout, is Tar more complicated. published on the sub- 
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: : : So : 
ject, surface perturbations have been depicted as being saw-tooth, ® continuous 


tatters ceseriped in the Serimger article, which are formed by wind and by 





SanvScid, anda randomly rough, in shape. In considering the extensive wave ! 


<ical variations, a logical refinement to the single-pulse model presented in 


mis work would be the extension of the concept to preduce a continuous wave 


ty 


omm (i.e., to represent the disturbed surface as a sinusoidal corrugation). 


In the case of Scrimger's bow-wave disturbance, the single-pulse concept 


ng is aa 


@ain 22 excended, this time to produce a wave train of Finite length travel 
ing sétweer source and receiver. Since, as noted earlier, the fluctuations are 
incuced mainiy by waves near the source and receiver, the question arises as to 


wnat portion of this elongated surface wave (the continuous wave or the extended 


bow-wave) should be considered in determining AI/I. A trial-and-error approach 
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| would allow one to find the areas where tne surface waves influence the 

| receiver, anc to estimate the effect of interference phenomenay sou, seme seca 
, prove tedicus. No attempt has been mace nere to solve this additional problem. 
} The sathematical formulation needed to develop the refinements outlined 
Mii@sov~ follows directly from the initial colLcests presented im Section ie oer 
puise can be rcslaced by a dipole, with the relative size of the pulse (i.e., 
wave neig¢gnt) determining the dipole strencth. Thus a decaying wave-form can be 


a secuence of dipoles of vaeryin 
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Sub cece wave does not have equal areas iying above and below the unperturbed 
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as indicated in the previous paragraphs. Instead, to an extent depend- 
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ing on the weve height, crests are relatively narrow and high whereas troughs 
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lativeiy broad and shallow. Tnus, to a first approximation, we could 
neglect tnose Cipoles lying below the true average surface, aliowing our model 
to consist of a series of "positive" dipoles, lying above a slightly lowered 
surface, and separated by a finite distance horizontally. 

B. Bottom Limitations -- Nowhere in the development of the mathematical 
model nor sudsequently in the results derived from it was mention made of the 
M@etccclons imposed Dy the presence of a lower boundary. In considering only 
the direcz ~ey and the once reflected surface ray (i.e., that ray produced by 
tne perturbztion dipole) a substantial simplification is injected into the 


meolem. According to theor ee lower boundary condition usually assumed is 
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| Duecan be a@monstrated that the chenge in exphLitude ard phase for a reflected 
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imumemee WAVS Will be given by the Rayleigh sezlection coefficients. ft canese 
® shown furcaer that this reflection coefficient is valid for spherical waves in 
i i 
I e - ~ ° = 2 - * vee ee Cex! ae ; 
Meee O61C" where the fractional variation in the space rate of change of the 
| Mol cuce secomes small 29 In the case under consideration the problem is com- | 
| + . - : ey > ) _ ake : = 
j polsacsed br ine fact that several sediments go into the make-up cf the bottom, | 
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EnG@ che tnic<iess associated with each layer varies along the source-receiver 
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‘| BS ef icitial attempt in refining the mod@l’to include lower boundamy 


erfecces, we: can consider the bottom to be pexvfectly lével and composed of a 
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mmeisee Section =). | 
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Simzle homogeneous material of infinite depth. If we thus assume the hottom to 
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be uniform, as well as "faster" and more dense than the water, the fact that the 
SsOurce and receiver lie well within @ wavelength of the bottom presents no probD- 
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lem. if the material beneath the source is considered rigid, it acts like an 


as ee a . e e 2 e 
He ite rigid baffle, and the energy radiated by the source is double that 
exzected in an unbounded medium. By the reciprocity theorem, the receiver now 
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WM critical dngle equals arcsin(c,/c,), where c, and ¢, are the propaga- 
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becomes doubly sensitive, so that the overall transmitter-receiver effect 
increases two-fold. However, since Tq. (3.1) is a ratio of intensities, the 


increase in the source strength cancels in the expression and does not affect 


ae von—a “5 ~ es Sa - 15  ~=a=— - Sa cen Sa sie Su = aaa, ” a tah} 
weit oO Cee ae eres ew @ a~a ete i Jo Neen a ret ha Se Ase aa Pe at 5+ 5 LLG Sue LS PiOwe SChnsc i en 


C. Velocity Gradient -- Throuzhour Secticas €1 and fi) iv WweseGaee 
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Cadac &4n Lsovelocity gradient existed trroughout the shallow water duct. 
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iheeect, oA1c Esquimalt Lagoon usually exhibits @ slignt positive gradicnc, sUseem 
Meecoxd ction of & positive grecient of scund speec, there is an Upward re7page 
ee oc; the “cys, anc because of this upward tending a larrer portion of the 
@e gy wilt Se vretlected from the bottom at angies of incidence sreater than 


tom reflection. 
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Bmome Cheveocteristics of the positive velocity gradient combine to reduce the 


Poss Of eftrsy from the shallow-water duct, thus producing, in general, better 
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MeeGemeso On characteristics. These effects ere, however, important only at 
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ane wise of 4900 ft/sec as the sound speed in Section Lil was appropriate 
Since Ref. ~ showed the profile to be isovelocity at this value, from a depth 
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EeemeoOxlnacely 3 ft down to tne bottom, In general, the procedure in céalmne 
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sc velocity gradients is to approximate the actual velocity pro- | 


File as a sepivc. of constant slopes, thus allowing for a cumulative solution of 


Gack. space nt. The partition of the velocity gradient in this case will do 


q 


actle to iitprove the results, since the scurce and receiver are bottom-mounted. 
fee ccc COI Wrd moveWnt of eny ray is quéth EZimited, end for all intents aad 


e607 


Mee oses tid WLocity profile listed in Ref. _ cah be considered isovelocity. 













D. The Dipole Assumption -- Assuming that tne acoustic effect of a surface 
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